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空间 通信 中 GMSK 信号 的 载波 恢复 与 时 钟 同步 
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摘 要 : 在 空间 通信 中 ， 过 大 的 多 普 勒 频 移 和 多 普 勒 变化 率 会 造成 相干 解 调 接收 机 无 法 进行 载波 恢复 与 时 钟 同步 。 
为 此 ， 针 对 高 斯 最 小 频 移 键 控 〈《GMSK) 信号 提出 一 种 可 靠 且 快速 克服 大 多 普 勒 频 移 的 改进 平方 环 方法 。 首 先 通过 
频率 估计 模块 获取 粗略 载波 频率 ， 其 次 应 用 改进 平方 环 结构 实现 GMSK 信号 的 载波 恢复 和 时 钟 同步 。 仿真 结果 显示 
在 输入 信 骂 比 Eb/No=5 dB 以 及 环 路 骂 声 带宽 为 50 kHz 的 条 件 下 ,载波 恢复 环 路 与 时 钟 同步 环 路 锁定 时 间 约 为 3.6 ms。 
在 低 信 骂 比 及 大 多 普 勒 频 移 的 条 件 下 , 该 方法 能 够 快速 且 稳 定 地 实现 GMSK 信号 在 空间 通信 中 的 载波 恢复 与 时 钟 同 
步 。 
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Carrier recovery and clock synchronization of gmsk signals in space communications 


Tang Zhiling, Zou Xin, Li Simini 
(Guangxi Key Laboratory of Wireless Broadband Communication & Signal Processing, Guilin University of Electronic 
Technology, Guilin Guangxi 541004, China) 


Abstract: Due to excessive Doppler shifts and Doppler rate in space communications, the coherent demodulation receiver 
fails to perform carrier recovery and clock synchronization. To this end, this paper proposed an improved square-loop 
method to overcome the large Doppler shifts reliably and quickly for the Gaussian Minimum Shift Keying (GMSK) signal. 
This method could obtain the coarse carrier estimation by the frequency estimation module, and then the carrier recovery 
and clock synchronization were realized by the improved square-loop. With loop noise bandwidth of 50 kHz and E/No=5 
dB, the simulation showed that the locked time of carrier-recovery loop and clock-synchronization loop was about 3.6 ms. 
In space communications , the method is proved can i1mplement carrier recovery and clock synchronization of GMSK 
signals quickly and stably under the conditions of low SNR and large Doppler shifts. 
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程 中 应 用 最 为 广泛 的 两 种 环 路 。 在 环 路 设计 中 环 路 捕获 带宽 
和 环 路 锁定 精度 是 相互 矛盾 的 。 尤 其 在 空间 通信 中 ， 接 收 信 
于 深 空 通信 任务 中 XX 波段 频谱 分 配 的 稀缺 性 , 深 空 通 号 中 往往 存在 较 大 的 多 普 勒 频 移 和 多 普 勒 变化 率 ， 如 火星 探 
信任 务 开始 转向 Ka 波段 进行 数据 通信 。 但 是 ， 在 深 空 通信 测 任务 中 多 普 勒 频 移 最 大 可 以 达到 90 kHz 左右 , 多 普 勒 变化 
任务 中 依然 需要 具备 X 波段 通信 系统 ， 用 作 应 急 通信 [5 。 率 最 大 达到 1000 Hz/st1。 如 果 想 要 提高 捕获 范围 大 小 ,高 
针对 频谱 资源 有 限 ， 数 据 传 输 速率 高 ， 传 输 距离 远 等 问题 ， 味 着 增加 环 路 噪声 带宽 , 在 低 信 噪 比 下 , 增加 环 路 噪声 带宽 ， 
高 斯 最 小 频 移 刍 控 (Gaussian minimum shift keying, GMSK) 会 造成 锁定 精度 下 降 。 因 此 ， 在 实际 应 用 中 需要 寻找 折 中 的 
调制 体制 以 其 恒 包 络 ， 频 谱 利 用 率 高 的 特点 被 国际 空间 数据 办 法 来 达到 工程 所 需求 的 性 能 。 与 此 同时 自动 频率 控制 环 路 ， 
系统 咨询 委员 会 (Consultative Committee for Space Data 同样 可 以 完成 某 些 调制 信号 的 载波 恢复 和 频 偏 估计 1, 通常 
Systems，CCSDS ) 列 为 深 空 通信 调制 方式 之 一 855。 在 2008 ”应 用 突 发 通信 ， 载 波 恢复 速度 很 快 ， 在 信 品 比 低 的 情况 下 时 
年 ， 欧 洲 航天 局 (European Space Agency, ESA) 将 GMSK 往往 锁定 精度 达 不 到 要 求 。 
调制 体制 首次 应 用 于 Herschel Planck 卫星 测试 任务 中 器。 在 实现 载波 恢复 与 频 偏 估计 的 工程 应 用 中 , 文献 [10, 11] 
GMSK 信和 号 的 解 调 方式 包括 非 相 干 解 调 与 相干 解 调 ， 传 中 提出 使 用 通用 数字 环 路 对 GMSK 调制 信号 进行 载波 恢复 ， 
统 的 GMSK 解 调 多 采用 非 相干 解 调 , 优点 是 接收 机 设计 复杂 但 是 由 于 鉴 相 器 具有 正弦 鉴 相 特性 ， 在 较 大 误差 信号 下 ， 鉴 
度 小 ， 但 是 解 调 性 能 比 相干 解 调 有 一 定 的 恶化 。 在 特定 领域 相 曲 线 并 非 呈 现 线性 特性 。 文 献 [12] 中 提出 基于 FPGA 线性 
如 空间 通信 方面 ， 由 于 非 相 干 解 调 存在 解 调 门限 的 限制 ， 只 鉴 相 全 数字 锁 相 环 ， 采 用 数字 计算 机 算法 (coordinate rotation 
能 通过 相干 解 调 方式 接收 GMSK 信号。 digital computer CORDIC) 实 现 线性 鉴 相 范围 的 扩展 。 而 针对 
在 GMSK 相干 解 调 接收 机 中 , 通常 采用 锁 相 环 结构 实现 空间 通信 中 出 现 的 大 多 普 勒 频 移 ;， 文献 [13] 提 出 一 种 基于 
相干 载波 的 提取 ， 其 中 平方 环 和 科斯 塔 斯 环 是 在 载波 恢复 过 FFT 算法 的 快速 且 精 准 的 载波 恢复 方法 ， 结 构 简 单 ， 估 计 精 
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信号 载波 恢复 与 时 钟 同步 的 FPGA 实现 ， 提 出 一 种 可 行 的 办 


法 。 
1 ”问题 分 析 


空间 飞行 器 ， 包 括 人 造 卫 星 、 载 人 飞船 、 航 天 飞机 等 ， 
在 轨 运 行 时 都 具有 很 快 的 运动 速度 。 因 此 ， 地 面 接收 的 信号 
中 必然 会 存在 着 多 普 勒 现象 。 多 普 勒 现象 包括 多 普 勒 频 移 以 
及 多 普 勒 变化 率 071， 两 者 的 数值 变化 大 小 ， 对 于 锁 相 环 参数 
设计 及 其 重要 。 
如 图 1 所 示 ， 以 近 地 轨 道 为 模板 ， 假 设 空间 飞行 器 绕 地 
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度 高 ， 通 过 频率 估计 加 快 锁 相 环 锁定 ， 文 献 [14,.15] 提 出 了 针 人 

对 QPSK 调制 信号 的 最 大 似 然 频 偏 估计 方法 ，FPGA 设计 较 lns 

为 复杂 ， 但 是 估计 精度 高 ， 文 献 [16] 提 出 在 FFT 频 偏 估计 上 pcos es 

进行 优化 算法 实现 精度 提升 ， 实 际 工程 应 用 较 难 。 Vrmsing (4) 
为 了 在 大 多 普 勒 频 移 与 低 信 品 比 条 件 下 实现 GMSK 信 0050 90o<g<l80s 

号 的 载波 恢复 与 时 钟 同步 ， 本 文 结合 前 人 的 工作 ， 提 出 将 平 Vt in0 

方 环 与 自动 频率 控制 环 路 相 结 合 的 改进 平方 环 方法 。 运 用 基 将 式 (3)(4) 代 入 式 (2) 中 可 以 得 到 

于 CORDIC 算法 的 FPGA 实现 完成 平方 环 的 线性 鉴 相 。 为 了 f 一 cosO 六 (5) 

实现 简易 而 精准 的 频 偏 估 计 , 采用 基于 FFT 算法 的 自动 控制 Vt nsinO c 

环 路 的 设计 ,辅助 平方 环 进行 载波 恢复 ,为 空间 通信 中 GMSK ”其 中 02<9<180° 。 


多 普 勒 变化 率 a 为 多 普 勒 频 移 在 单位 时 间 1 内 的 变化 
值 ， 多 普 勒 频 移 f 对 时 间 + 进行 求 导 可 以 得 到 多 普 勒 变化 率 
公式 为 


= 


a,= 
“~ 


—vV,, (rn, sin 0—(r2 + )sinO+rn,) 大 


V2 十 用 2 —2n,n, sin 0) C 


(0) 


其 中 0"<O<180。 。 
如 图 1 所 示 ， 当 飞行 器 处 于 弧 线 SS" 范围 内 才 可 以 进行 
正常 的 信号 发 送 与 接收 。 即 


T 


球 进 行 近似 圆周 运动 且 地 面 接收 机 与 飞行 器 轨道 处 于 同一 平 
面 。 根 据 飞行 器 运动 轨迹 图 分 析 地 面 接收 机 与 空间 飞行 器 通 
言 中 存在 的 多 普 勒 现象 。 


图 1 飞行 器 运动 轨迹 平面 图 
Fig.1 Aircraft motion trajectory plan 

以 地 心 0 为 原点 建立 坐标 系 , 其 中 R 点 为 地 面 接收 机 位 

置 ，S 点 为 空间 飞行 器 当前 位 置 ， 地 球 半径 为 上 ， 轨 道 高 度 

为 n，S 点 处 的 线 速 度 为 w ， 假 设 飞行 器 从 B 点 运行 到 S 点 


处 ， 转 动 角度 为 2。 4 是 RS 与 08 的 夹 角 。 由 牛顿 万 有 引力 
公式 可 得 
d0 GM 
a 3 
d G+ (1) 
| Vs = Oh 


中 :=i+h 为 飞行 器 轨道 半径 ，% 是 飞行 器 的 角速度 ，G 
万 有 引力 常量 , M 为 中 心 星体 的 质量 , 此 处 是 地 球 的 质量 。 
根据 多 普 勒 频 移 计算 公式 求 得 多 普 勒 频 移 fo 为 
fa = 

其 中 : vw 为 飞行 器 线 速度 v, 在 向 量 RS 的 速度 分 量 ，/ 为 无 线 
电信 号 传输 频率 ，c 为 光 传 播 速率 。 
速度 分 解法 可 得 

Ie —v,sing 0°<0<90° 

w= 90° <0<180° G3) 


并 六 


站 


时 根据 三 角 函 数 定理 可 得 


可 


B=arcsin(r,/r,) 
Ls 0) 
基于 以 上 公式 的 推导 ，X 波段 无 线 电信 号 传输 频率 /为 
8.4 GHz, 通过 改变 飞行 器 的 高 度 , 对 多 普 勒 现象 的 特性 进行 
分 析 。 
如 图 2 所 示 , 飞行 器 高 度 越 低 , 多 普 勒 频 移 最 大 值 越 大 。 
高 度 为 300 KM 时 , 多 普 勒 频 移 最 大 值 可 达到 200 kHz 左右 ， 
在 0=90* 时 多 普 勒 频 移 达到 最 小 值 0 Hz。 


200 T T 


一 1500KM| 
一 2000KM 


= 

© 

= 
T 


普 勒 频 移 (kHz) 


多 
© 
已 


40 60 80 100 120 14 
度 8(") 
到 2 ”接收 信号 中 的 多 普 勒 频 移 
Fig.2 Doppler shifts of receive signals 
如 图 3 所 示 ， 飞 行 器 高 度 越 低 ， 变 化 率 最 大 值 也 越 大 。 
高 度 为 300 KM 多 普 勒 变化 率 最 大 值 可 达 3 kHz 左右 ， 且 在 
9=90°? 时 多 普 蔓 变化 率 取得 最 大 值 。 
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度 8(") 
图 3 接收 信号 中 的 多 普 勒 变化 率 


Fig.3 ”Doppler rate of receive signals 
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2 GMSK 载波 恢复 与 时 钟 同步 方法 


GMSK 信号 是 将 码 元 信号 预先 经 过 一 个 窗 带 高 斯 滤波 
器 ， 然 后 再 进行 MSK 调制 ， 它 同样 可 以 解释 为 一 种 调制 指 
数 为 0.5 的 连续 相位 的 频 移 调制 (CPFSK) 081。 鉴 于 在 无 线 
通信 中 采用 的 是 抑制 载波 调制 方式 ， 虽 然 接收 信和 号 
载波 的 频率 分 量 ， 但 是 信号 含有 载 频 信息 ， 只 要 经 过 非 线性 
变换 即 可 产生 载波 的 倍 频 分 量 。 
2.1 GMSK 信号 分 析 
根据 MSK 调制 的 特点 ， 每 个 码 元 周期 相位 变化 为 
tmt/127 〈 允 是 码 元 周期 )， 那 么 GMSK 信号 中 每 个 码 元 周期 


的 频率 变化 范围 为 
= 上 -1147 
全 =f.+1/4, (8) 
其 中 大 为 载波 频率 。 
假设 输入 信号 为 
x(t)=s(t)+n(7) 
-msn| ost 如 |+n() (9) 
其 中 : ml(?) 是 信号 幅度 失真 ，n(?)=n.(D)+ jn(1)，n.( 站 与 (7) 
是 完全 独立 的 高 斯 白 噪 声 ， 为 载波 角 频 率 。 


GMSK 信号 平方 后 经 带 通 滤波 器 得 


y0=F3mO rn Om -d+ 
cos[2@.t+20,(7)]+ 
[m(7)—n, (7)]n.(t) sin[2wz+26 (7)] 
@ (1)=+xt/27, 。 
式 (10) 知 平方 后 信号 存在 离散 频率 2fs 和 2f., GMSK 
盲 号 平方 频谱 如 图 4 所 示 。 
GMSK 信 号 二 次 方 的 频谱 


(10) 


H.C ; 
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妈 5 ”改进 平方 环 结构 图 


Fig.5 Modified square loop structure diagram 


2.2.1 全 数字 锁 相 环 模型 分 析 


6 是 二 阶 模拟 锁 相 环 相位 模型 图 ， 鉴 相 器 鉴 相 增益 为 
Ks， 坏 路 滤波 器 S 域 传递 函数 为 F(9)， 压 控 振荡 器 $ 域 传递 


函数 为 V(s) ， 鉴 相 器 


输出 误差 电压 V(s), 经 环 路 滤波 器 滤波 


得 控制 电压 V.(s) ， 控 制 输出 相位 8(s) 靠拢 输入 相位 4(s) 。 


模拟 锁 相 环 相位 模型 


在 线性 鉴 相 区 间 的 S 域 变换 公式 为 


_V(s)_1lrstl 


MO Ts 0 

Vs)= 人 -全 (13) 
H(s) @(s) KF(s)V(s) 
A(s) TCFCJYO 

250s + (14) 


其 中 : 


S72 +260,s + ,7 


K, 为 压 控 振荡 器 (VCO) 增 益 ，7 和 亏 是 时 间 常 数 ， 


8 
@,=NKiK, /nn 为 环 路 然 频率 ， <=zo,/2 为 阻尼 因子 。 
6 - 
HH(s) 为 模拟 锁 相 环 在 S 域 的 传递 函数 。 
是 a 9.(s) K, Vi (SF la 
世 
哩 27 一 
0 AR 二 | 加 2 Se 
od 图 6 二 阶 模拟 锁 相 环 相位 模型 
-3 2 1 加 0 1 2 3 Fig.6 Second-order analog phase-locked loop phase model 
频率 (Hz) x10 为 了 对 模拟 锁 相 环 数字 化 ， 可 以 通过 双 线 性 变换 方法 Pal 
图 4 ”GMSK 信号 平方 频谱 实现 模拟 系统 与 数字 系统 之 间 的 转换 。 变 换 公式 为 
Fig.4 Spectrum of GMSK signal square _21-2 
载波 频率 与 离散 频率 之 间 的 关系 为 1 
二 ts (15) 
7 (11) 一 
J =2fu -2f 2 
获取 2/ 和 2 两 个 离散 频率 信号 , 混 频 后 经 低 通 滤波 就 。 其 中 民 是 数字 采样 周期 。 
可 以 得 到 时 钟 频率 为 信号， 经 脉冲 形成 后 ， 即 可 提取 时 钟 将 式 (15) 分 别 代入 式 (12)(13) 得 
言 号 。 1 
) = C+ C 一 -一 1 
2.2 改进 平方 环 分 析 与 设计 ST 0 
数字 平方 环 是 数字 锁 相 环 的 一 种 ， 包 括 鉴 相 器 、 环 路 滤 N= Re (17) 
波 器 和 数字 控制 振荡 NCO)。 平 方 环 的 设计 需要 分 析 数 字 Lk 
锁 相 环 和 模拟 锁 相 环 来 确定 环 路 参数 ， 文 献 [19] 对 基于 二 阶 其 中 : QG=(2m+7)/2r ， C=T /no 
模拟 锁 相 环 的 数字 锁 相 环 设计 作 了 详细 的 分 析 。 本 文 提 出 的 数字 锁 相 环 的 相位 模型 如 图 7 所 示 。 则 数字 锁 相 环 z 域 
改进 平方 环 结 构 如 图 5 所 示 。 和 专递 函数 可 以 表示 为 
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KCzt+(KC, — KC )z? 
1+(KC —2)z'+(KC, — KC +1)z2 


其 中 KK=KsK。， 为 数字 环 路 总 增益 。 
A 0.(z) K Va) Eta (mca) 


图 7 二 阶 数字 锁 相 环 相 位 模型 

Fig.7 Second-order digital phase-locked loop phase model 

式 (18) 是 根据 数字 锁 相 环 相位 模型 推导 出 来 的 z 域 系统 
传递 函数 ， 另 外 还 可 以 将 模拟 锁 相 环 系 统 传递 函数 ， 采 用 双 
线性 变换 方法 变换 成 数字 域 的 系统 函数 , 将 式 (15) 代 入 式 (14) 
可 以 得 到 


H(z)=(460,T,. +2(@,T) |+2(@T) z+ 
2(@T) zth(o,T) -460,T |z?)x 
|(@T -4s0,T |z*)x(4+4s0,T.+(2T) |+ 


[2(@,7) -8 +[4-4s0,T +(@,T) |] (19) 


通过 比较 式 (18) 和 (19) 可 以 得 出 
4(@T) +8é0T, 1 
"4+4e0,T +(@,T)Y K 
AoT) 1 
~ 4+4é@T, +(oT) K 
2.2.2 改进 平方 环 参数 计算 
式 (18)(20) 可 知 ， 数 字 锁 相 环 系统 传递 函数 中 ,需要 计 
算 鉴 相 增益 K, 、 数 字 控 制 振荡 器 增益 到 、 环 路 总 增益 KK、 
阻尼 系数 8 、 环 路 自然 角 频 率 w 、 二 阶 环 路 滤波 器 积分 系数 
G 与 比例 系数 C,。 
首先 ， 文 献 [21] 中 提出 当 阻 尼 系 数 了 =0.707 以 及 @,T. <<1 
时 ， 工 程 应 用 效果 最 好 。 为 采样 周期 。 
1) 环 路 自然 角 频 率 w 
数字 锁 相 环 捕获 过 程 包括 频率 捕获 与 相位 捕获 ， 改 进 平 
方 环 采用 自动 控制 环 路 对 大 频 偏 进行 估计 ， 所 以 这 里 讨论 相 
位 捕获 ， 快 捕 带 被 称 为 保证 环 路 只 有 相位 捕获 过 程 的 最 大 固 
有 频率 值 。 故 快 捕 带 为 


(20) 


Aw, =2é0, (21) 
为 满足 环 路 稳定 环 路 自然 角 频率 w 必须 满足 

MW (22) 
针对 频率 斜 逢 信号， 频率 变 化 率 为 R 时， 二 阶 锁 相 环 路 


维持 环 路 稳定 ， ww 满足 公式 : 
四 >VR (23) 


同时 ， 还 要 从 环 路 噪声 性 能 上 对 @ 的 进行 分 析 ， 在 空间 
通信 中 接收 信号 存在 高 斯 白 噪声 的 影响 ， 可 以 得 到 在 输入 信 


号 中 存在 高 斯 白 噪声 条 件 下 相位 方差 可 以 表示 为 

和 ma 
其 中 : (S1N) 为 输入 信 噪 比 ，B 为 输入 信号 带宽 ， 有 为 环 路 
呆 声 带宽 。 
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环 路 噪声 带宽 越 小 ,相位 方差 o? 也 就 越 小 ， 稳 态 波动 也 
就 越 小 , 对 于 理想 二 阶 环 路 滤波 器 环 路 噪声 带宽 B 可 以 表示 
为 


OO， 
B=ge(+46) (25) 


环 路 信 噪 比 (85/N), 其 定义 为 环 路 输入 信号 功率 与 单 边 

带 到 的 噪声 功率 之 比 ， 表 示 环 路 对 噪声 抑制 能 力 ， 通 常 当 
(S/N),>6dB 时 环 路 才 可 以 正常 锁定 ，(3/N), 可 以 表示 为 

[国光 芭 09 


则 由 以 上 分 析 可 以 得 到 


国民 
”NA3(+422) 
根据 式 (22)(23) 和 (27) 可 以 确定 环 路 自然 角 频 率 w 的 范 


(27) 


Am 有 < 人 总 ) 4éB, 

22 "WN)3(1+46) 
wT <l1 (28) 
四 >VR 


2) 鉴 相 增 益 K， 

鉴 相 器 是 一 个 相位 比较 装置 ， 用 来 检测 输入 相位 信号 与 
反馈 信号 之 间 的 相位 差 ,反正 切 鉴 相 器 工作 原理 如 图 8 所 示 。 

输入 信号 请 = | 误差 信号 


图 8 反正 切 鉴 相 原理 图 
Fig.8 Arctangent Phase discrimination schematic 
假设 鉴 相 器 输入 信号 为 
u(t)=U,cos(@t+0 (1)) (29) 
中 : w 为 输入 信号 角 频 率 ，4.(7) 为 输入 信号 初始 相位 。 
数字 控制 振荡 器 输出 本 振 信 号 可 以 表示 为 
有 号。 cos(w,t +0,(1)) 
wo lt)=U, sin(@t+0,(1t)) 
其 中 : w 为 本 振 信 号 角 频 率 ， 4(0) 为 本 振 信 号 初始 相位 ， 环 
路 锁定 时 四 =w 。 
输入 信号 与 本 振 信号 经 过 混 频 ， 
2 倍 频 信号 后 得 


并 


(30) 


过 低 通 滤波 器 滤 除 


[E: 


UU, 


(0=2 eola ()] 
(31) 
zuo (1)= sin[8, (1)] 


其 中 : 4.(7)=8(7)-8(?) 。 则 反正 切 鉴 相 器 输出 误差 信号 表达 式 


us (t)=arctan (wa (1)/ uzo (t)) 


=K,0.(1) (32) 
即 反 正切 鉴 相 增益 Ks =1。 
3) 数 字 控 制 振荡 器 增益 K。 
数字 控制 振荡 器 它 的 作用 就 是 产生 正 交 的 正弦 余弦 样 
本 ， 工 作 原理 如 图 9 所 示 。 
根据 锁 相 环 工作 原理 ， 环 路 滤波 器 输出 的 控制 电压 用 以 
修正 数字 控制 振荡 器 的 输出 频率 ， 则 NCO 输出 频率 为 


f= 坟 (Wh +AW) G3) 


201904.00009v1 


国 
国 


XIV 


china 


唐 智 灵 ， 等 : 
其 中 ， 为 采样 频率 ，N 为 相位 累加 字 字 长 ，W 为 初始 频率 
字 ，A 为 控制 频率 字 ， 数 字 控 制 振荡 器 相位 增益 为 


2 
K, = ET (34) 


中 Ti 为 相位 累加 字 更 新 周期 工程 应 用 中 取 Ti =87 。 


并 


累加 器 


图 9 数字 控制 振荡 器 原理 图 
Fig.9 Numerically controlled oscillator schematic 
如 图 5 所 示 的 改进 平方 环 结构 ， 输 入 信号 由 于 位 宽 的 变 
化 会 增加 额外 的 增益 ， 即 环 路 总 增益 由 环 路 滤波 器 输出 位 宽 
8 所 影响 ， 最 终 环 路 滤波 器 增益 K 可 以 重新 表示 为 
K-20 ) KK, = ET 247) (35) 
空间 通信 中 ， 信 号 处 理 都 是 将 高 频 信 号 搬移 到 中 频 进行 
数据 处 理 ， 本 文中 数字 信号 采样 率 人 A =50 MHz， 载 波 频率 
大 =10 MHz， 信 息 速率 妨 =1 MHz， 则 GMSK 输入 信号 带宽 
B; =500 kHz 。 
根据 图 2 和 3 的 多 普 勒 现象 分 析 ， 多 普 勒 频 移 最 大 200 
kHz， 多 普 勒 变化 率 最 大 3 kHz。 根 据 载 波 估计 频 偏 精度 ， 快 
捕 带 Aw. 最 大 为 3 kHz， 频 率 变化 率 R 最 大 为 3 kHz, 输入 信 
品 比 乌 /No。=5dB， 本 文 已 知 阻尼 因子 了 =0.707， 鉴 相 增 益 
Ks=1， 为 了 提升 锁 相 环 精度 频率 字 位 宽 N=32， 环 路 滤波 器 
位 宽 轧 =29， 根 据 式 (28) 和 (35) 兼 顾 锁 定时 间 与 稳 态 相差 则 可 
以 设 定 @ 为 


@, =2xx5x103 (rad/s) (36) 
根据 式 (35)(26) 及 (18) (20) 可 以 计算 改进 平方 环 z 域 传 


0.0267z- 一 0.0263z 一 
1 一 1.9733z- +0.9737z”™ 


根据 式 (37) 可 以 计算 改进 平方 环 的 极点 为 0.986730.0132i ， 
极点 位 于 圆 内 则 证 明了 该 环 路 是 稳定 的 。 
2.3 CORDIC 鉴 相 的 FPGA 实现 

CORDIC 算法 可 以 将 三 角 函 数 计算 转换 为 加 减 和 移 位 操 
作 ， 这 非常 适合 于 在 硬件 电路 上 实现 ， 该 算法 不 仅 可 以 节约 
FPGA 的 硬件 资源 ， 还 能 大 大 提高 计算 速度 [ 久 。CORDIC 算 
法 坐标 旋转 图 如 图 10 所 示 。 


H(z)= (37) 


了 


(x yin) 


图 10 ”坐标 旋转 图 
Fig. 10 ”Coordinate rotation diagram 


旋转 公式 为 
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Xu =XcOsO—y,sing 
i (38) 
yi =N SinO+y;cosO 
通过 去 除 因数 cosb 可 以 得 到 伪 旋 转 公式 为 
X= ytan0=x—y27 
| p -i (39) 
yn = y+xtan0= y+x27 


HH tan0=27i。 
同时 伪 旋 转 公式 可 以 表示 为 
Xn =X — y27d;? 
yi = td (40) 
Zin = 2 —d;arctan27 

其 中 : z 是 迭代 角度 的 累加 值 ， 每 次 迭代 旋转 的 角度 为 
arctan 2 ， 大 小 可 以 通过 查 表 法 获取 , 取 值 如 表 1 所 示 。d; 等 
于 sgn〈yY)， 决 定 每 次 旋转 的 方向 。 
表 1 旋转 角度 值 


Table 1 Rotation angle value 
i tan0=27 0/(°) 
1 1.00000000 45.00000000 
之 0.50000000 26.56505118 
3 0.25000000 14.03624347 
4 0.12500000 7.12501635 
5 0.06250000 3.57633437 
6 0.03125000 1.78991061 
7 0.01562500 0.89517371 
8 0.00781250 0.44761417 
9 0.00390625 0.22381050 
CORDIC 鉴 相 曲线 如 图 11 所 示 。 
41 
8 
五 
全 
以 
-4 ! ! ! | 
-10 -5 0 5 10 
输入 相位 (rad) 
了 11 CORDIC 算法 鉴 相 曲线 
Fig. 11 CORDIC algorithm phase discrimination curve 


2.4 载波 频率 估计 的 FPGA 实现 

锁 相 环 路 的 捕获 时 间 与 初始 频 差 有 着 直接 的 关系 ， 初 始 
频 差 越 大 ， 捕 获 时 间 越 长 ， 受 限于 锁 相 环 路 的 性 能 要 求 ， 捕 
获 带 宽 较 小 ， 则 在 频 偏 较 大 的 情况 下 无 法 完成 本 地 载波 与 术 
位 的 同步 ， 继 而 无 法 进行 相干 解 调 。 通 过 对 输入 信号 进行 分 
析 确 定 粗 略 载波 频率 ， 控 制 初始 频 差 在 一 个 较 小 的 范围 内 ， 
便于 锁 相 环 路 进行 快速 捕获 。 
频率 估计 模块 的 原理 框图 如 图 12 所 示 。 


JU 


过 


输入 有 频率 字 
频率 /in 再 上 加 N 点 谱 线 向 频率 2 in, 频率 字 | 输出 、 
再 FFT 估计 


机 


tt SSH 
tt 


洪 


12 频率 估计 原理 图 
Fig. 12 Frequency estimation schematic 
频率 估计 模块 的 核心 是 图 12 中 的 N 点 的 FFT 运算 ( 采 
用 中核 )。 通过 对 NN 个 输入 信号 点 进行 FFT 运算 后 ,对 运算 
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结果 进行 求 模 并 查找 最 大 幅 值 的 位 置 ， 即 谱 线 点 ， 通 过 谱 线 ”输入 信 噪 比 为 为 /mm =20 dB 和 如 /N,。=5 dB， 进 行 ModelSim 


点 的 位 置 确定 估计 值 ， 最 后 转换 成 频率 字 输 出 。 仿真 如 图 15 和 16 所 示 。 
本 方法 是 基于 FFT 算法 实现 频率 估计 ， 结 构 简 单 ， 对 输 i 


入 信 品 比 为 轧 /N =5dB， 频 移 在 -300 kHz~300 kHz 范围 内 的 
GMSK 信和 号 的 估计 精度 进行 了 仿真 如 图 13 所 示 , 结 果 显 示 能 
够 满足 多 普 勒 频 移 估 计 要 求 。 | 

x10 , 
3 上 EbIN0=5dB| 


* .Pll_vig_tst/din |11001000 


A 


/ 
WwW 


人 


“poy e0002795520|8... ,880279552 


fb 
T 


二 
T 
\ 


.oll_vig_tst/b_clk IStO 


四 | 
15 ”载波 恢复 与 时 钟 同步 ( /No =20dB) 
Fig. 15 Carrier recovery and clock synchronization ( E, /No =20db) 


-3 Ea ee 
- -2 -1 0 1 2 3 


实际 载波 频 移 值 (Hz) x105 
图 13 FFT 载波 估计 仿真 图 


Fig. 13 FFT carrier estimation simulation diagram 
SU 人 \ 十 、 
3 ”实验 与 讨论 


文中 选用 Altera 公司 的 Cyclone IV 系列 的 
EP4CE75F23C8 芯片 ， 采 用 QuartusII 和 ModelSim 以 及 
MATLAB 进行 联合 仿真 测试 传统 平方 环 与 改进 平方 环 在 大 


™ 


pa 


载波 频 移 估计 值 (Hz) 


' 
Dm 


"Pll_vlg_tst/din |00101110 


多 普 勒 频 移 ， 以 及 不 同 信 品 比 条 件 下 的 工作 性 能 。 rm 
测试 输入 初始 条 件 : 调制 方式 为 GMSK， 采 样 率 为 50 Bow ooo00ps 
MHz， 载 波 频率 为 10 MHz， 初 始 频 移 为 250 kHz， 多 普 勒 变 图 16 ”载波 恢复 与 时 钟 同步 ( /No =5dB) 
化 率 为 3 kHz/s， 信 息 速率 1 Mbps。 Fig. 16 Catrrier recovery and clock synchronization( E, /No =5db) 
3.1 传统 平方 环 结构 15 中 ,载波 恢复 环 路 的 环 路 滤波 器 输出 frequency_df _c 


传统 平方 环 中 只 包含 鉴 相 器 、 环 路 滤波 器 、 数 字 控 制 振 。” 和 时 钟 同步 环 路 的 环 路 滤波 器 输出 frequency_df_s 同时 进入 
荡 器 (NCO) 三 个 主要 模块 ， 其 中 鉴 相 器 一 般 采 用 的 正弦 鉴 ”” 收 仇 。 说 明 在 输入 信 噪 比 为 By/N =20 dB 且 频 差 处 于 较 小 范 
相 器 实现 近似 线性 鉴 相 。 内 的 情形 下 , 传统 平方 环 可 以 完成 GMSK 信号 的 载波 恢复 
实验 1 传统 平方 环 结构 进行 GMSK 调制 信号 的 载波 恢 “ 与 时 钟 同步 。 载波 频 率 估计 模块 在 1.2 ms 左右 获取 载波 频率 

复 与 时 钟 同步 ， 输 入 信 品 比 为 鸟 /N=20 dB, 进行 ModelSim ”估计 的 频率 字 f_c ， 将 频 差 控制 在 2197 Hz 左右 ， 环 路 自然 
仿真 如 图 14 所 示 。 角 频 率 为 50 kHz 的 传统 平方 环 在 5.5 ms 左右 完成 载波 环 路 

. 与 时 钟 环 路 的 锁定 。 

16 中 ， 信 噪 比 为 为 /Nm =5 dB 的 情况 下 ， 传 统 平方 环 
ce 载波 恢复 环 路 的 环 路 滤波 器 输出 fequency_df _c 和 时 钟 同步 
“Pl va terst | | 环 路 的 环 路 滤波 器 输出 frequency_df _s 无 法 收敛 。 环 路 未 能 实 
-LVO_ydout 0000111000 EE 现 载波 恢复 与 时 钟 同步 的 功能 ， 说 明 在 低 信 品 比 条 件 下 传统 


.frequency df.c -1278977 Te 平方 环 无 法 维持 其 环 路 的 稳定 性 。 
Pantano 3.2 改进 平方 环 结构 
| 本 文 提出 改进 平方 环 的 结构 是 在 平方 环 基础 上 增加 了 载 


波 估 计 模 块 应 对 大 多 普 勒 频 移 影响 ， 减 小 初始 频 差 大 小 ， 同 


SS 


fr 


.pl vig tstdk | 


-LVg_tststarf_c|857714482 时 利用 CORDIC 算法 的 FPGA 实现 完成 线性 鉴 相 , 减少 鉴 相 
E00 2 误差 带 来 的 影响 ， 用 以 提升 锁 相 环 的 稳定 性 。 
图 14 载波 恢复 与 时 钟 同步 ( /No =20dB) 实验 3 ”改进 平方 环 结构 ， 输 入 信 噪 比 为 为 /No =20 dB， 
Fig. 14 Carrier recovery and clock synchronization ( E, /No =20db) 进行 ModelSim 仿真 如 图 17 所 示 。 
图 14 中 ,载波 恢复 环 路 的 环 路 滤波 器 输出 frequency_df -ce 将 图 17 中 环 路 滤波 器 输出 、 载 波 输出 dout 和 时 钟 信号 
和 时 钟 同步 环 路 的 环 路 滤波 器 输出 frequency_df_s 的 输出 结  b_clk 等 数据 进行 MATLAB 分 析 ， 结 果 如 图 18 和 19 所 示 。 


果 显 示 未 收 化 ， 说 明 环 路 处 于 失 锁 状态 ， 证 明 在 大 多 普 勒 频 。 ”经 过 1.2 ms 左右 的 频 偏 估计 ， 将 初始 频 偏 控制 在 2197 kHz 


de 


移 和 快速 多 普 勒 变化 率 的 情况 下 ， 环 路 自然 角 频 率 为 50kHz ”左右 ， 环 路 自然 角 频 率 为 50 kHz 的 改进 平方 环 在 2 ms 左右 
时 ， 由 于 环 路 捕获 带宽 小 于 固定 频 差 ， 传 统 平方 环 无 法 实现 ”完成 相位 锁定 ， 并 输出 频率 为 10.25 MHz 的 相干 载波 ， 图 20 
GMSK 信和 号 的 载波 恢复 与 时 钟 同 步 。 可 得 提取 的 时 钟 信号 周期 为 1 MHz。 说 明 在 输入 信 噪 比 为 

实验 2 ”在 传统 平方 环 的 基础 上 加 入 载波 估计 模块 ， 加/No=20 dB 及 接收 信号 中 存在 大 多 普 勒 频 移 和 快速 多 普 勒 
的 是 克服 大 多 普 勒 频 移 的 影响 ,将 初始 频 差 控 制 在 较 小 范围 ， ”变化 率 的 情况 下 ， 短 时 间 内 实现 载波 恢复 与 时 钟 同 步 。 
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图 17 载波 恢复 与 时 钟 同步 ( EE, /No =20 dB) 
Fig. 17 Carrier recovery and clock synchronization( E, /No =20 dB) 
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图 18 ”载波 环 路 仿真 数据 图 ( /No =20 dB) 
Fig. 18 Carrier loop simulation data map( E, /No =20 dB) 
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图 19 ”载波 输出 仿真 数据 图 ( 5B, /No =20 dB) 
Fig. 19 Carrier output simulation data map( &, /No =20 dB) 
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图 20 ”时钟 同步 环 路 仿真 数据 图 ( /No =20 dB) 


Fig. 20 Clock synchronization loop simulation data 
diagram( E, / No =20 dB) 
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实验 4 改进 平方 环 结构 ,输入 信 噪 比 为 为 /MX =5 dB, 进 
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行 ModelSim 仿真 如 图 21 所 示 。 
将 图 21 中 环 路 滤波 器 输出 、 载 波 输出 dout、 时 钟 信号 


b_clk 等 数据 进行 MATLAB 分 析 ， 结 果 如 图 22 和 23 所 示 。 
载波 频 偏 估计 同样 将 频 差 控 制 在 2197 Hz 左右 ， 与 输入 信 噪 


比 为 轧 /No=5dB 的 情况 下 相 比 ， 频 偏 估计 时 间 延 长 至 2.2 ms 
左右 ， 环 路 锁定 时 间 延 长 至 3.6 ms 左右 ， 输 出 频率 为 10.25 
MHz 的 相干 载波 ， 图 24 可 得 时 钟 信 号 输出 周期 为 1 MHz。 


此 验证 在 输入 信 品 比 为 马 /N,=5dB 及 接收 信号 存在 大 多 


定 工 作 。 


普 勒 频 移 和 多 普 勒 变化 率 的 情况 下 ， 改 进 平 方 环 依然 可 以 稳 
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图 21 载波 恢复 与 时 钟 同步 ( /No =5 dB) 


Fig. 21 Carrier recovery and clock synchronization( E/No =5 dB) 
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图 22 载波 环 路 仿真 数据 图 ( /No =5 dB) 


Fig. 22 Carrier output simulation data map( E, /No =5 dB) 
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图 23 载波 输出 仿真 数据 图 ( E/No =5 dB) 


Fig. 


23 Carrier output simulation data map( E, /No =5 dB) 


综合 以 上 实验 可 知 传统 平方 环 无 法 在 大 多 普 勒 频 移 、 快 


荔 变 化 率 以 及 低 信 品 比 条 件 下 保证 环 路 正常 工作 的 稳 


定性 ， 


站 | 


本 文 提 出 的 改进 平方 环 方案 能 够 在 大 多 普 勒 频 移 及 


低 信 品 比 的 条 件 下 完成 GMSK 信号 的 载波 恢复 与 时 钟 同步 ， 


与 此 同时 


[在 环 路 的 锁定 速度 与 稳定 性 上 也 有 


定 的 提高 。 
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